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Evolution  of  the  Ion  Velocity  Distribution  in  the 
Near  Field  of  the  BHT-200-X3  Hall  Thruster 
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Edwards  AFB,  CA  93524 
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Abstract 

This  work  presents  an  analysis  of  near  plume  velocity  distributions  of  the  Busek  BHT-200-X3  200  W  laboratory  Hall  thruster 
derived  from  laser-induced  fluorescence  measurements  of  the  5 d[4] 7/2- 6p [3] 5/2  xenon  ion  excited  state  transition.  The  evolution  of  the 
axial  and  radial  velocity  distributions  is  presented  from  the  exit  plane  to  approximately  0.7  exit  diameters  downstream.  It  is  shown  that 
the  distributions  evolve  significantly  between  the  exit  channel  center  and  the  inner  portion  of  the  channel  where  there  are  high  frac¬ 
tions  of  low  velocity  ions.  It  is  believed  that  this  may  be  due  to  the  slight  inward  focus  of  the  xenon  ion  propellant  colliding  on  the 
surface,  recombining,  and  subsequently  being  re-ionized.  Low  axial  velocity  ions  appear  downstream  of  the  thruster  exit  plane.  While 
their  origin  is  not  entirely  certain,  the  most  likely  explanation  appears  to  be  cross  annular  flow.  Strong  evidence  of  charge  exchange  is 
also  evident  producing  high  velocity  wings  on  the  primary  ion  velocity  peak  reminiscent  of  ion  energy  analyzer  studies  of  other  thrus¬ 
ters.  A  survey  of  the  axial  and  radial  velocity  distributions  along  the  exposed  exterior  boron  nitride  nose  cone  covering  the  central 
magnetic  pole  is  performed  to  determine  the  behavior  of  the  ions  impacting  this  critical  surface.  The  VDF  is  ideal  for  comparison  with 
numerical  simulations  since  it  allows  for  direct  comparison  of  the  fundamental  ion  acceleration.  External  VDF  data  can  even  provide 
limited  insight  as  to  the  internal  ion  formation  and  acceleration  processes. 


Introduction 

The  goal  of  this  study  is  to  characterize  the  xenon  ion  ve¬ 
locity  distribution  function  (VDF)  of  200  W  Busek  BHT-200- 
X3  xenon  Hall  thruster  in  the  region  between  the  exit  and  cath¬ 
ode  planes.  The  VDF  is  ideal  for  comparison  with  numerical 
simulations  since  it  allows  for  direct  comparison  of  the  funda¬ 
mental  ion  acceleration.  External  VDF  data  can  even  provide 
limited  insight  as  to  the  internal  ion  formation  and  acceleration 
processes. 

Many  previous  attempts  to  compare  bulk  ion  velocities 
from  experiments  with  results  of  numerical  simulations  have 
encountered  difficulties.  For  example,  simulation  efforts  usu¬ 
ally  report  mean  velocities.  Experimental  methods  such  as 
laser-induced  fluorescence  (LIF)  are  often  noise  limited  and  are 
better  suited  to  determining  the  most  probable  velocity.  In  the 
skewed,  non-symmetric  velocity  distributions  common  to  Hall 
thrusters,  the  most  probable  (e.g.  peak  signal)  and  statistical 
mean  velocities  differ.  Depending  on  the  shape  of  the  distribu¬ 
tion,  these  two  bulk  measures  of  velocity  may  vary  signifi¬ 
cantly.  The  advantage  of  using  a  single  bulk  velocity  measure  is 
that  it  allows  for  the  succinct  presentation  of  large  data  sets. 
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However,  both  the  modeler  and  experimentalist  are  discarding 
the  majority  of  the  generated  data.  Comparison  of  the  VDF 
determined  from  each  method  will  allow  for  minimal  ambigu¬ 
ity.  Furthermore,  the  annular  geometry  of  the  Hall  thruster  dis¬ 
charge  /  acceleration  chamber  produces  a  complex  flow  field 
where  multiple  velocity  populations  interact.  The  comparison 
of  VDF  data  is  these  cases  is  vital  to  increasing  understanding 
of  these  flows. 

Xenon  Laser  Induced  Fluorescence 

LIF  is  a  convenient  diagnostic  for  the  investigation  of  the 
VDF  as  it  does  not  perturb  the  plasma.  The  LIF  signal  is  a  con¬ 
volution  of  the  VDF,  transition  line  shape,  and  laser  beam  fre¬ 
quency  profile.  Determination  of  the  VDF  from  LIF  data  only 
requires  the  deconvolution  the  transition  line  shape  and  laser 
beam  profile  from  the  raw  LIF  signal. 

For  the  results  reported  here,  the  5 d[4]  7/2- 6p  [3]  5/2  electronic 
transition  of  Xe  II  at  834.7  nm  is  probed.  The  isotopic  and 
nuclear-spin  effects  contributing  to  the  hyperfme  structure  of 
the  5J[4]7/2-6/7[3]5/2  xenon  ion  transition  produce  a  total  of  19 
isotopic  and  spin  split  components.  The  hyperfme  splitting  con¬ 
stants  which  characterize  the  variations  in  state  energies  are 
only  known  for  a  limited  set  of  energy  levels.  Unfortunately, 
the  834.7  nm  xenon  ion  transition  only  has  confirmed  data  on 
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the  nuclear  spin  splitting  constants  of  the  6p[3]s/2  upper  state 
[1-4].  Manzella  and  others  have  previously  used  the  5cf[4]7/2- 
6p[3]5/2  xenon  ion  transition  at  834.7  nm  to  make  velocity 
measurements  in  a  Hall  thruster  plume  [1,5-7].  A  convenient 
feature  of  this  transition  is  the  presence  of  a  relatively  strong 
line  originating  from  the  same  upper  state  (6S[2]3/2-6p[3]s/2 
transition  at  541.9  nm  [8]  which  allows  for  non-resonant  fluo¬ 
rescence  collection. 

Ion  velocity  is  determined  by  measurement  of  the  Doppler 
shift  of  the  absorbing  ions.  If  an  absorber  has  a  velocity  com¬ 
ponent  u  along  the  axis  of  the  laser  beam,  it  will  absorb  the 
light  at  a  frequency  shifted  from  that  of  stationary  absorbers. 
The  magnitude  of  this  frequency  shift  dvu  depends  on  the  speed 
u  along  the  laser  beam  axis 


dV^ 


[1] 


where  c  is  the  speed  of  light.  The  Doppler  shift  of  a  species’ 
fluorescence  profile  away  from  the  line  center  Vn  of  stationary 
absorbers  is  in  proportion  to  u  [9]. 

Experimental  Apparatus 

The  measurements  in  this  effort  were  performed  in  Cham¬ 
ber  6  of  the  Air  Force  Research  Laboratory  (AFRL)  Electric 
Propulsion  Laboratory  at  Edwards  AFB,  CA.  Chamber  6  is  a 
non-magnetic  stainless  steel  chamber  with  a  1.8  m  diameter  and 
3  m  length.  It  has  a  measured  pumping  speed  of  32,000  1/s  on 
xenon.  Pumping  is  provided  by  four  single  stage  cryogenic 
panels  (single  stage  cold  heads  at  ~25  K)  and  one  50  cm  two 
stage  cryogenic  pump  (<12  K).  Chamber  pressure  during  thrus¬ 
ter  operation  is  approximately  7xl0"4  Pa,  corrected  for  xenon. 


Fig.  L  Side  view  diagram  of  thruster  within  AFRL 
Chamber  6.  Also  shown  are  the  translation  stages 
and  the  laser  probe  beams.  Note  that  the  fluores¬ 
cence  collection  and  external  optics  are  not  shown. 


The  Hall  thruster  used  in  this  study  is  the  Busek  Company 
BHT-200-X3  200  W  Hall  thruster  which  has  been  described 
elsewhere  [10].  Table  1  shows  the  nominal  operating  conditions 
for  the  BHT-200-X3  thruster  during  this  test.  During  thruster 
operation,  the  parameters  shown  in  Table  1  are  monitored  and 
recorded  at  a  1  Hz  data  rate. 


Table  1.  Nominal  Thruster 

Operating  Conditions 

Anode  Flow 

840  (Ig/s  (Xe) 

Cathode  Flow 

98  pg/s  (Xe) 

Anode  Potential 

250  V 

Anode  Current 

0.82  A 

Keeper  Current 

0.5  A 

Magnet  Current 

1.0  A 

Heater  Current 

3.0  A 

Figure  1  shows  a  side-view  diagram  of  the  Hall  thruster 
mounted  within  the  vacuum  chamber.  The  thruster  is  mounted 
on  a  three  axes  orthogonal  computer  controlled  translation  sys¬ 
tem.  Figure  1  also  shows  the  two  orthogonal  LIF  probe  beams 
and  windows  through  which  the  beams  enter  the  chamber.  Fig¬ 
ure  2  shows  a  top  view  of  the  laser  optical  train,  collection  op¬ 
tics,  through  which  the  beams  end  one  leg  of  the  external  probe 
optics. 

The  laser  used  is  a  New  Focus  Vortex  tunable  diode  laser. 
It  is  capable  of  tuning  approximately  ±50  GHz  about  a  center 
wavelength  of  834.7  nm.  The  10  mW  beam  is  passed  through  a 
Faraday  isolator  to  eliminate  feedback  to  the  laser.  The  laser 
beam  then  passes  through  several  beam  pick-offs  until  it 
reaches  a  50-50  beam  splitter  (BS)  where  it  is  split  into  two 
beams  of  equal  power.  The  first  beam,  the  axial  probe  beam 
shown  in  Figs.  1  and  2,  is  focused  by  a  lens  and  enters  the  vac¬ 
uum  chamber  through  a  window.  A  second  probe  beam,  shown 
in  Fig.  1  only,  is  directed  from  the  optical  bench  via  a  periscope 
apparatus  so  that  it  enters  the  chamber  from  above  the  thruster 
and  probes  the  velocity  perpendicular  to  the  first  probe  beam. 
Each  probe  beam  is  chopped  at  a  unique  frequency  by  choppers 
Ch2  (2.8  kHz)  and  Ch3  (2.0  kHz)  for  simultaneous  phase  sensi¬ 
tive  detection  of  the  two  fluorescence  signals. 

The  two  wedge  beam  pick-offs  (/BS)  shown  in  Fig.  2  pro¬ 
vide  portions  of  the  beam  for  diagnostic  purposes.  The  first 
beam  pick-off  directs  a  beam  to  a  photodiode  detector  (Dl) 
used  to  provide  constant  power  feedback  to  the  laser.  The  sec¬ 
ond  beam  is  divided  into  two  equal  components  by  a  50-50 
cube  beam  splitter.  The  first  component  is  directed  to  an  Bur¬ 
leigh  WA-1500  wavemeter  used  to  monitor  absolute  wave¬ 
length.  The  second  component  is  sent  through  a  chopper  Chi 
(1.3  kHz)  and  through  a  low  pressure  xenon  hollow  cathode 
discharge  lamp.  The  lamp  provides  a  stationary  absorption  ref¬ 
erence  for  the  determination  of  the  Doppler  shift  dvu.  Unfortu¬ 
nately,  there  is  no  detectable  population  of  the  ionic  xenon 
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collection  optics.  Note  that  the  radial  probe  beam  peri¬ 
scope  and  focusing  optics  are  not  shown. 

5 of [4] 7/2  state.  However,  there  is  a  nearby  (measured  to  be  17.6 
GHz  distant  using  this  apparatus)  neutral  xenon  6s'  [1/2] i-6p ' [3/ 
2]2  transition  at  834.68  nm  [11-12].  The  second  pick-off  sends  a 
beam  to  a  300  MHz  free  spectral  range  Fabry-Perot  etalon  (F- 
P)  that  provides  high  resolution  frequency  monitoring  of  the 
wavelength  interval  swept  during  a  laser  scan. 

The  fluorescence  collection  optics  are  also  shown  in  Fig.  2. 
The  fluorescence  is  collected  by  a  75  mm  diameter,  300  mm 
focal  length  lens  within  the  chamber.  The  collimated  fluores¬ 
cence  signal  is  directed  through  a  window  in  the  chamber  side 
wall  to  a  similar  lens  that  focuses  the  collected  fluorescence 
onto  the  entrance  slit  of  125  mm  focal  length  monochrometer 
with  a  photomultiplier  tube  (PMT).  Due  to  the  1:1  magnifica¬ 
tion  of  the  collection  optics,  the  spatial  resolution  of  the  meas¬ 
urements  is  determined  by  the  geometry  of  the  entrance  slit  1 
mm  width  and  1.5  mm  height  as  well  as  the  sub-mm  diameter 
of  the  probe  beams. 

Figure  3  shows  the  near  field  geometry  of  the  Busek  BHT- 
200  Hall  thruster.  The  locations  of  the  protruding  central  mag¬ 
netic  pole  (nose  cone)  and  edges  of  the  acceleration  channel  are 
indicated  as  is  the  position  of  the  cathode  exit.  The  cartesian 
coordinate  system  and  origin  used  in  these  measurements  are 
also  shown  in  Fig.  3.  The  coordinate  system  orientation  is  also 
referenced  in  Figs.  1  and  2.  The  origin  is  0.5  mm  beyond  the  tip 
of  the  nose  cone  due  to  the  repeatability  with  which  this  posi¬ 
tion  may  be  located.  All  locations  are  referenced  to  these  coor¬ 
dinates.  Measurements  presented  in  this  work  will  be  limited  to 
the  Y-Z  plane.  Therefore  in  the  anode  discharge  exit  region,  the 
Y  axis  velocity  measurements  correspond  to  the  radial  veloci¬ 
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Fig.  3.  Near  field  dimensions  of  the  BHT-200-X3 
Busek  Co.  Hall  thruster  with  origin  of  the  coordinate 
system  and  positions  of  critical  dimensions  and  loca¬ 
tions  noted. 

ties,  and  the  Z  axis  measurements  to  the  axial  velocities. 

Velocity  Distribution  Extraction 

The  convolution  of  two  functions  a(t )  and  z(t)  is  defined  as 

ja(z)z{t-z)dz  =  a(t)*z{t)  =  u(t) 

J  P] 

where  *  denotes  convolution  and  u(t )  is  the  resulting  convolved 
function.  In  the  specific  case  of  interest  here,  u(t)  is  the  raw  LIF 
profile,  z(t)  is  the  VDF,  and  a{t)  is  the  convolution  of  the  transi¬ 
tion  line  shape  and  any  other  profile  broadening  mechanisms 
such  as  laser,  or  instrument  broadening  and  is  usually  denoted 
in  the  literature  as  the  point  spread  function  (PSF)  [13]. 

For  analytical  functions,  solution  of  Eqn.  2  is  trivial.  Dis¬ 
crete  data  with  a  reasonably  high  signal  to  noise  ratio  (SNR) 
may  employ  a  number  of  numerical  techniques  with  mixed 
results.  The  introduction  of  the  noise  inherent  to  the  measure¬ 
ments  will  usually  create  a  mathematically  ill-posed  problem 
that  lacks  a  unique  solution.  Noise  tends  to  produce  instabilities 
in  the  numerical  deconvolution  process  which  propagate 
through  the  solution  and  contaminate  the  results.  Deconvolu¬ 
tion  then  becomes  a  tradeoff  between  solution  smoothness  and 
numerical  accuracy.  Subsequently,  the  deconvolution  of  noisy 
data  may  become  subjective.  Iterative  methods  attempt  to  miti¬ 
gate  this  issue  by  incorporating  physical  constraints  in  a  rigor¬ 
ous  methodology  with  noise  analysis  to  determine  the  extent  of 
deconvolution  possible.  However,  the  effect  of  noise  in  real 
data  during  deconvolution  limits  the  application  of  all  deconvo¬ 
lution  algorithms  [13-14]. 
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For  discrete  data,  convolution  may  be  represented  by 

Az-u  p] 

where  u  is  the  convolved  vector  with  noise  inherent  to  the  data 
collection,  A  is  the  PSF  in  matrix  format,  and  z  is  the  desired 
vector,  in  this  case  the  VDF.  The  deconvolution  algorithm  used 
in  this  work  is  a  zeroth  order  Tikhonov’s  regularization  method 
[14-15].  Tikhonov’s  regularization  method  seeks  to  solve  Eqn. 
3  through  use  of  an  artificial  smoothing  function  to  cure  the  ill- 
conditioned  nature  of  the  deconvolution.  It  does  this  by  mini¬ 
mizing  an  artificial,  constrained  smoothing  function  M\z,u) 
defined  to  be 

Ma (z,u)  =  \\Az - uf  -  «||z||2 


where  a  is  the  regularization  parameter.  The  minimization  is 
performed  by  setting  the  derivative  of  Ma(z,u)  to  zero  which 
yields 

(ATA  +  al)z  =  ATu  ^ 

where  I  is  the  identity  matrix.  Equation  5  may  be  solved  for  z. 
In  this  work,  Eqn.  5  is  solved  for  z  using  LU  decomposition. 

The  precept  of  Tikhonov  regulation  is  that  a  unique,  stable 
solution  to  the  deconvolution  may  be  achieved  by  minimizing 
the  constrained  smoothing  function  Ma(z,u).  Fundamentally, 
this  will  generate  a  stable,  single  valued  solution  only  so  long 
as  the  uncertainty  of  the  iterated  solution  does  not  exceed  the 
uncertainty  (noise)  of  the  original  data.  Within  these  con¬ 
straints,  the  algorithm  is  relatively  insensitive  to  noise  [15]. 

The  choice  of  a  affects  the  deconvolution  process.  Too 
high  of  a  value  will  overly  smooth  the  result  and  too  low  of  a 
value  will  allow  the  noise  produce  nonphysical  results.  Happily, 
regularization  provides  a  mechanism  to  determine  the  degree  to 
which  deconvolution  can  proceed  [16].  The  L-curve  is  a  log-log 
plot  of  the  norm  of  the  residual  versus  that  of  the  corresponding 
norm  of  the  regularized  solution.  As  the  regularization  parame¬ 
ter  is  varied,  the  L-curve  provides  a  convenient  graphical  tool 


Fig.  4.  Hollow  cathode  discharge  lamp  UF  pro¬ 
file  of  xenon  ion  5cf[4]7/2-6p[3]5/2  transition  and 
gaussian  curve  fit  used  as  the  PSF. 
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for  determining  the  trade  off  between  regularized  solution  and 
fit  to  the  data.  The  vertex  of  the  L-curve  denotes  the  best  trade 
off  between  these  two  sources  of  uncertainty  by  balancing  the 
magnitude  of  the  solution  and  the  solution  smoothness. 

The  main  advantage  of  the  L-curve  method  is  that  it  allows 
for  a  quantitative  methodology  to  determine  the  optimal  regu¬ 
larization  parameter  a.  If  an  estimate  of  the  noise  inherent  to 
the  data  is  available,  it  is  possible  to  determine  similar  criterion 
without  generating  the  multiple  deconvolutions  required  to  con¬ 
struct  the  L-curve.  However,  experience  has  shown  that  noise 
estimates  are  not  sufficiently  accurate  and  that  use  of  the  L- 
curve  in  determining  the  optimum  regularization  parameter 
eliminates  many  of  the  uncertainties  in  the  noise  estimation. 
Most  importantly,  the  L-curve  removes  a  potential  source  of 
experimental  bias  from  the  data  analysis. 

For  this  effort,  the  PSF  of  the  LIF  data  was  measured  using 
the  apparatus  shown  in  Fig.  2  with  several  minor  modifications. 
Fluorescence  was  collected  from  the  low  pressure,  pass  through 
hollow  cathode  lamp  (<100  Pa)  using  the  125  mm  mono¬ 
chrometer  with  a  single  300  mm  focal  length  lens.  A  typical 
LIF  line  profile  is  shown  in  Fig.  4  along  with  a  Gaussian  curve 
fit.  Although  the  transition  line  shape  is  not  strictly  symmetri¬ 
cal,  a  Gaussian  profile  with  a  full  width  at  half  maximum 
(FWHM)  of  600  MHz  fits  the  data  adequately.  This  PSF  ac¬ 
counts  for  both  instrument  and  line  shape  effects. 

Figure  5  shows  a  sample  axial  LIF  profile  and  resulting 
deconvolved  VDF  from  the  Hall  thruster  channel  exit  plane 
where  the  raw  LIF  profile  is  broadest.  The  data  is  an  average  of 
5  raw  LIF  traces  of  1000  points  filtered  with  5  passes  of  a  25 
point  second  order  Savitzky-Golay  filter  interpolated  to  512 
points  for  the  deconvolution.  The  data  has  an  estimated  SNR  of 
50.  The  PSF  is  substantially  narrower  than  the  broad  LIF  trace 
with  an  equivalent  FWHM  of  500  m/s.  As  a  result,  deconvolu¬ 
tion  does  not  substantially  change  this  raw  LIF  trace.  Figure  6 
shows  a  radial  distribution  located  28  mm  downstream  from  the 
axial  VDF  in  Fig.  5.  Here  deconvolution  narrows  the  peak,  in¬ 
creasing  the  peak  height  by  approximately  20%  (Note  that  both 
raw  and  deconvolved  data  are  area  normalized).  The  deconvo¬ 
lution  of  this  LIF  profile  results  in  a  greater  change  than  that 
seen  in  Fig.  5  due  to  the  narrower  peak  in  the  raw  data.  Broader 
features  such  as  that  between  2-5  km/s  are  not  affected  by  the 
deconvolution.  It  should  also  be  noted  that  the  deconvolution 
routine  typically  produces  oscillatory  noise  near  the  trace  end¬ 
points  as  seen  in  Fig.  5  and  cropped  out  in  Fig.  6. 

Deconvolution  is  not  strictly  required  to  estimate  xenon  ion 
VDFs  from  the  raw  LIF  data  in  this  particular  plasma  dis¬ 
charge.  Not  performing  the  deconvolution  introduces  an  uncer¬ 
tainty  conservatively  estimated  to  be  less  than  20%.  The  major¬ 
ity  of  the  following  data  will  consist  of  data  that  has  been  de¬ 
convolved.  However,  some  data  with  low  SNR  with  broad  dis¬ 
tributions,  especially  those  near  the  exit  plane,  will  not  be  de¬ 
convolved  to  minimize  deconvolution  induced  noise.  All  such 
data  will  be  clearly  identified. 
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Fig.  5.  Exit  plane  axial  LIF  profile  (U  vector),  and  re¬ 
sulting  deconvolved  VDF  (Z  vector).  Also  shown  is  the  L- 
curve  computed  with  the  location  of  regularization  pa¬ 
rameter  used.  Note  the  small  difference  between  VDF  and 
LIF  profile 


Fig.  6.  Radial  LIF  profile  (U  vector)  28  mm  downstream 
of  exit  plane,  and  resulting  deconvolved  VDF  (Z  vector). 
Also  shown  is  the  L -curve  computed.  Note  the  sharpen¬ 
ing  of  the  peak  due  to  deconvolution  is  significant  due  to 
the  relatively  narrow  distribution  relative  to  the  P5F. 

Results  and  Discussion 

Understanding  and  characterizing  the  near  plume  region  is 
important  in  understanding  the  interaction  of  the  ion  plume 
with  the  environmental  background  gas  as  well  as  a  starting 
point  for  the  modeling  of  the  plume  expansion  to  determine  the 
affects  of  thruster  placement  on  thruster/spacecraft  interaction. 
It  has  been  shown  previously  that  the  xenon  ions  emitted  from  a 
Hall  thruster  undergo  approximately  50  eV  of  additional  accel¬ 
eration  between  the  thrusters  physical  exit  plane  and  cathode 
plane  [17].  The  ionic  flow  evolves  in  the  near  plume  region  and 
the  thruster  can  not  be  fully  characterized  without  studying  this 
region  in  detail.  In  this  specific  thruster  with  its  protruding  cen¬ 
tral  magnetic  pole,  the  erosion  of  the  boron  nitride  nose  cone  is 
a  life  limiter.  Knowledge  of  the  ionic  VDF  in  this  region  would 
aid  in  the  prediction  of  thruster  life  and  provide  guidance  con¬ 
structing  longer  lived  Hall  thrusters. 

Channel  Center  VDF  Evolution 

Figures  7  and  8  show  the  evolution  of  the  axial  and  radial 
VDF  with  varying  distance  from  the  exit  plane  (Z  =  -7  mm) 
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through  the  cathode  plane  (Z  =  +13  mm)  to  a  distance  28  mm 
from  the  exit  plane  (Z  -  +21  mm).  All  data  in  Figs.  7  and  8  is 
from  the  center  of  the  acceleration  channel  (Y  =  12  mm).  In 
each  figure,  the  VDF  data  is  area  normalized  and  all  8  traces 
are  provided  with  identical  scaling.  In  both  axial  and  radial 
data,  the  VDF  retains  a  fundamental  form  as  the  ions  travel 
through  the  near  plume  region.  However,  there  are  a  number  of 
important  changes  that  provide  insight  of  the  processes  occur¬ 
ring  to  the  ions  as  they  travel  through  the  near  plume. 

The  axial  VDF  typically  consists  of  a  narrow  primary  peak 
overlaying  a  broader  distribution  of  generally  low  velocity 
classes.  The  main  peak  is  accelerated  between  the  exit  plane 
and  nose  cone  tip  (Z  =  0)  by  approximately  2.1  km/s,  from 
15,100  to  17,200  m/s.  The  peak  continues  to  accelerate  beyond 
the  nose  cone  tip,  but  not  as  significantly.  The  central  peak  also 
narrows  during  the  external  acceleration  as  each  ion  gains  ap¬ 
proximately  47  eV  from  the  exterior  electric  field.  With  the 
corresponding  velocity  increment  dependent  on  initial  velocity, 
the  central  peak  sharpens.  The  most  interesting  portion  of  the 
axial  VDF  evolution  along  the  channel  centerline  is  the  appear¬ 
ance  of  populations  of  low  speed  ions.  As  expected,  the  broad 
low  energy  portion  of  the  ion  distribution  is  depleted  during 
between  Z  =  -7  mm  and  Z  =  +1  as  the  ions  are  accelerated. 
However  at  Z  =  +5  mm,  a  new  class  of  low  axial  velocity  ions 
appears.  There  continues  to  be  populations  of  low  axial  velocity 
ions  appearing  beyond  Z  >  5  mm. 

The  radial  VDF  data  presented  in  Fig.  8  shows  a  generally 
sharp,  nearly  symmetrical  distribution.  Similar  to  the  axial  data, 
the  radial  VDF  peaks  narrow  with  increasing  distance  along  the 
Z  axis.  This  can  be  explained  by  the  flow  divergence,  fixed 
sample  volume,  and  increasing  distance  from  the  source.  Closer 
examination  of  the  data  reveals  that  there  is  a  slight  asymmetry 
in  the  data.  A  positive  velocity  tail  clearly  appears  just  beyond 
the  nose  cone  (Z  =  1  mm).  The  main  radial  peak  continues  to 
narrow  until  Z  =  17  mm,  despite  the  positive  velocity  tail.  By  Z 
=  21  mm,  this  tail  feature  feature  extends  to  +10  km/s,  ap¬ 
proximately  80  eV,  and  has  slightly  depressed  the  VDF  peak 
height. 

Figure  9  shows  the  measured  axial  VDF  at  the  exit  plane 
as  well  as  a  calculated  energy  distribution  and  the  integrals  of 
each.  The  axial  VDF  main  peak  is  located  at  15.1  km/s  (155 
eV)  with  a  width  of  approximately  4  km/s.  The  mean  axial  ve¬ 
locity  and  kinetic  energy  are  calculated  to  be  12.4  km/s  and  139 
eV.  Approximately  75%  of  the  flow  is  in  the  main  peak  and  it 
contains  approximately  90%  of  the  energy.  Interestingly,  the 
small  feature  at  20.1  km/s  has  a  peak  energy  of  290  eV  which  is 
nearly  twice  that  of  the  mean  and  peak  energies  (139-155  eV). 
Although  the  SNR  is  poor,  the  values  are  consistent  with  charge 
exchange  ions.  The  radial  exit  plane  VDF  data  in  Fig.  10  show 
a  narrow  distribution  with  a  symmetrical  distribution  within  a 
slightly  negative  mean  velocity  -570  m/s  (1.3  eV). 

Figures  11  and  12  show  the  corresponding  axial/radial 
VDF  data  at  Z  =  +5  mm.  The  main  axial  peak  has  accelerated 
to  16.9  km/s  (195  eV).  The  mean  axial  velocity  has  increased  to 
14.0  km/s  and  the  mean  kinetic  energy  is  180  eV.  A  low  veloc¬ 
ity  axial  peak  (-2  km/s)  first  appears  at  this  location  accounting 
for  13%  of  the  flow  but  less  than  2%  of  the  axial  energy  com- 
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Fig.  7.  Axial  VDF  functions  for  acceleration 
channel  centerline  (Y  =  12  mm)  from  the  exit 
plane  (Z  =  -7  mm)  to  28  mm  beyond  (Z  =  21 
mm).  Note:  VDFs  are  scaled  and  normalized 
to  unity  area. 


Fig.  8.  Radial  VDF  functions  for  acceleration 
channel  centerline  (Y  =  12  mm)  from  the  exit 
plane  (Z  =  -7  mm)  to  28  mm  beyond  (Z  =  21 
mm).  Note:  VDFs  are  scaled  and  normalized  to 
unity  area. 
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Fig.  9.  Centerline  axial  VDF,  inset  energy  distribution , 
and  relative  summations  (green)  at  the  thruster  exit 
plane  (Z  =  -1  mm).  Mean  axial  velocity  and  kinetic 
energy  are  12,400  m/s  and  139  eV,  respectively. 


Fig.  10.  Centerline  radial  VDF,  inset  energy  distribu¬ 
tion,  and  relative  summations  (green)  at  the  thruster 
exit  plane  (Z  =  -1  mm).  Mean  radial  velocity  and  ki¬ 
netic  energy  are  -570  m/s  and  -1.3  eV,  respectively. 

ponent.  The  main  axial  peak  still  contains  approximately  75% 
of  the  flow  and  90%  of  the  energy  deposited  in  the  flow.  Inter¬ 
estingly  enough,  the  axial  VDF  shows,  a  small  (-2%)  popula¬ 
tion  of  ions  with  velocities  greater  than  19.17  km/s  (250  eV). 
This  is  greater  than  the  applied  250  V  discharge  potential  and 
may  be  evidence  of  charge  exchange  collisions  in  the  plume. 
The  radial  distribution  in  Fig.  12  shows  noticeable  asymmetry 
with  a  positive  velocity  tail.  This  is  especially  evident  in  the 
energy  distribution  which  shows  20%  of  the  radial  energy  to  be 
in  this  asymmetric  feature. 

Figures  13  and  14  show  the  corresponding  axial/radial 
VDF  data  at  Z  =  +2 1  mm.  The  main  axial  peak  has  accelerated 
to  17.2  km/s  (202  eV).  The  main  axial  peak  now  only  contains 
65%  of  the  flow  and  75%  of  the  flow  energy.  The  mean  axial 
velocity  is  now  14.8  km/s  with  a  mean  axial  kinetic  energy  of 
186  eV.  Now  approximately  6%  of  the  ions  have  energies 
greater  than  250  eV.  This  increase  can  be  seen  more  globally  in 
Fig.  7  as  the  high  velocity  tail  of  the  main  axial  feature  grows 
with  Z.  The  radial  VDF  in  Figure  14  shows  a  larger  fraction  of 
the  radial  flow  (-15%)  is  now  in  the  high  velocity  tail  which 
clearly  extends  to  80  eV  and  that  it  now  accounts  for  approxi¬ 
mately  65%  of  the  radial  energy. 

The  asymmetry  beyond  the  nose  cone  is  believed  to  be  due 
to  mixing  with  the  flow  from  the  opposite  portion  of  the  thrus¬ 
ter  annulus.  The  behavior  of  the  radial  data  is  consistent  with 
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Fig.  11.  Centerline  axial  VDF,  inset  energy  distribu¬ 
tion,  and  relative  summations  (green)  at  Z  =  5  mm. 
Mean  axial  velocity  and  kinetic  energy  are  -13,985  m/ 
s  and  180  eV,  respectively.  Note  the  appearance  of  a 
low  velocity  ion  population. 


Fig.  12.  Centerline  radial  VDF,  inset  energy  distribu¬ 
tion,  and  relative  summations  (green)  at  Z  =  5  mm. 
Mean  radial  velocity  and  kinetic  energy  are  400  m/s 
m/s  and  6.5  eV,  respectively. 

the  apparently  low  and  mid  range  velocity  populations  in  the 
axial  VDF  data,  especially  if  collisions  between  the  various 
populations  create  low  velocity  populations.  It  also  appears  to 
be  consistent  with  peak  velocity  data  from  further  in  the  plume 
taken  in  a  previous  study  which  show  a  distinct  second  peak 
appearing  further  downstream  which  is  clear  evidence  of  cross 
annular  flow  [17].  The  axial  VDF  traces  in  Fig.  7  do  show  one 
feature  of  the  axial  VDF  data  which  can  not  be  explained  by 
cross  annular  flow.  Beginning  just  beyond  the  exit  plane,  the 
high  velocity  tail  originating  on  the  main  peak  becomes  notice¬ 
able  and  progressively  more  prominent  with  distance  from  the 
thruster  exit  plane.  This  appears  to  be  evidence  of  charge  ex¬ 
change  collisions  between  singly  and  doubly  charged  xenon 
ions  as  previously  reported  by  King  in  the  SPT-100  thruster 
using  parallel  plate  energy  analyzers  in  the  distant  plume  [18]. 
The  appearance  of  the  very  low  speed  ions  beyond  Z  =  1  mm 
may  also  be  related  to  charge  exchange,  but  of  collisions  of 
singly  charged  ions  with  neutrals.  However,  charge  exchange  as 
the  primary  source  low  energy  ions  remains  speculative  at  this 
time  and  it  is  believed  that  annular  cross  flow  with  significant 
collisionallity  is  the  more  likely  explanation. 

VDF  Evolution  Near  Thruster  Nose  Cone 

From  the  acceleration  channel  centerline  data  presented  in 
Figs.  7-14,  it  is  apparent  that  the  ion  beam  is  focused  slightly 
inward.  Since  Hall  thruster  lifetime  is  limited  by  the  erosion  of 
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Fig.  13.  Centerline  axial  VDF,  inset  energy  distribu¬ 
tion,  and  relative  summations  (green)  at  Z  =  21  mm. 
Mean  axial  velocity  and  kinetic  energy  are  14,750  m/s 
and  186  eV,  respectively. 


Fig.  14.  Centerline  radial  VDF,  inset  energy  distribu¬ 
tion,  and  relative  summations  (green)  at  Z  =  21  mm. 
Mean  radial  velocity  and  kinetic  energy  are  1,230  m/s 
and  19  eV,  respectively. 

the  insulator  channel  walls  and  for  this  thruster  also  that  of  the 
boron  nitride  nose  cone  covering  the  protruding  central  mag¬ 
netic  pole,  it  is  important  to  understand  interaction  of  surfaces 
with  the  plasma.  After  an  ion  strikes  a  wall,  it  is  typically  as¬ 
sumed  that  full  accommodation  occurs.  The  ion  recombines 
with  a  wall  electron  to  form  a  neutral  and  subsequently  leaves 
the  wall  fully  thermalized  without  a  preferential  direction. 
These  resultant  neutrals  may  then  interact  with  the  plasma  in  a 
variety  of  ways.  These  include  creating  regions  of  higher  than 
expected  density  with  correspondingly  higher  electron  conduc¬ 
tion  due  to  increased  numbers  of  collision  partners.  These  re¬ 
combined  ions  may  also  be  re-ionized  and  accelerated  by  the 
external  electric  fields  to  energies  of  40-50  eV  (7. 5-8. 5  km/s). 
Due  to  the  divergence  of  Hall  thruster’s  external  electric  fields, 
these  low  energy  ions  will  likely  be  accelerated  with  signifi¬ 
cantly  more  divergence  than  the  main  flow.  Understanding  this 
portion  of  the  flow  is  therefore  of  great  importance  to  the 
spacecraft  integrator. 

Figure  15  shows  a  survey  of  the  axial  and  radial  ionic  VDF 
data  along  the  protruding  central  magnetic  nose  cone  of  the 
thruster.  This  data  has  not  been  deconvolved  due  to  the  broad 
LIF  profiles  and  the  low  SNR.  However,  all  data  in  Fig.  15  has 
been  area  normalized  and  all  graphs  contain  a  common  scaling. 
Figure  15  also  shows  the  relative  size  and  placement  of  meas¬ 
urement  regions  as  well  as  the  calculated  mean  axial  and  radial 
velocities.  These  traces  are  from  a  single  series  of  LIF  traces 
and  contain  no  averaging,  but  have  been  filtered  to  remove  high 
frequency  noise.  The  flow  begins  at  the  exit  plane  (Z  =  -7  mm) 
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with  an  axial  VDF  with  a  high  proportion  of  low  velocity  ions. 
These  are  likely  a  result  of  the  central  nose  cone  geometry  that 
places  this  measurement  region  in  the  lee  of  the  interior  portion 
of  the  thruster  nose  cone.  The  radial  VDF  at  this  position  shows 
a  divergent  flow  toward  the  insulator  with  a  peak  at  -4  km/s.  In 
Fig.  15,  the  general  trend  of  the  axial  VDF  data  that  the  higher 
velocity  peak  grows  as  the  slow  moving  ions  are  externally 
accelerated.  This  trend  ends  at  Y,  Z  =  6,  -4  mm.  Beginning  here 
and  extending  to  the  nose  cone  tip,  the  axial  VDF  is  increas¬ 
ingly  dominated  by  the  low  velocity  portion  of  the  VDF.  The 
mean  axial  velocity  drops  steadily  from  this  point  onward  and 
the  distributions  broaden  significantly.  At  the  nose  cone  tip,  the 
axial  VDF  is  an  extremely  broad  feature  with  a  mean  velocity 
of  1  km/s.  It  is  interesting  to  note  that  after  some  period  of  run¬ 
ning  the  thruster  used  in  this  study,  an  indentation  formed  on 
the  nose  cone  tip.  The  broad  distribution  at  the  nose  cone  tip  (Y, 
Z  =  0,  0  mm)  in  Fig.  15  shows  that  nearly  half  the  ions  in  this 
region  are  moving  toward  the  thruster  and  are  subsequently 
sputtering  material  that  would  ordinarily  be  considered  to  be  in 
the  shadow  of  the  nose  cone. 

The  radial  distributions  in  Fig.  15  show  a  steadily  broaden¬ 
ing  distribution  along  the  nose  cone  profile.  The  radial  VDF 
peak  velocity  also  steadily  decreases  along  nose  cone;  starting 
at  -4  km/s  at  the  exit  plane  and  approaching  -8  km/s  just  prior 
to  the  tip.  At  the  nose  cone  tip,  the  radial  distribution  is  nearly 
20  km/s  wide  and  likely  contains  ions  from  both  sides  of  the 
annulus  which  contribute  to  the  width  of  the  distribution. 

Interestingly  at  the  exit  plane  (Z  =  -7  mm),  there  is  a  weak 
radial  VDF  peak  at  +4  km/s  that  appears  to  contain  approxi¬ 
mately  1 0%  of  the  apparent  population.  This  may  be  due  to  the 
reflection  of  the  laser  off  the  nose  cone  producing  a  spurious 
signal.  This  effect  has  been  used  before  by  Cedolin,  et  al.,  who 
used  reflection  off  a  Hall  thruster  anode  to  eliminate  the  need 
for  a  stationary  reference  source  [19].  However,  we  can  not 
completely  rule  out  that  the  signal  represents  a  distinct  popula¬ 
tion  of  ions.  This  illustrates  one  of  the  difficulties  inherent  to 
laser  diagnostics  for  plasma/surface  interactions.  For  the  analy¬ 
sis,  the  possibly  spurious  signal  is  retained  for  calculation  of 
mean  velocities  since  the  possibility  of  ions  flowing  in  away 
from  the  nose  cone  surface  can  not  be  completely  ruled  out, 
especially  further  along  the  nose  cone.  At  the  nose  cone  tip,  it 
appears  that  there  are  ions  from  both  sides  of  the  annulus  which 
produce  the  extremely  broad  axial  and  radial  distributions. 
There  are  some  similarities  at  Z  =  -6  mm  where  there  is  also  an 
apparent  positive  velocity  flow.  It  may  be  that  the  small  posi¬ 
tive  velocity  portion  of  the  radial  VDF  is  due  to  leakage  around 
the  nose  cone  of  plasma  from  the  opposite  portion  of  the  annu¬ 
lus.  It  is  also  important  to  consider  the  dimensions  of  the  sam¬ 
ple  volumes  which  may  explain  these  results. 

Exit  Plane  Profiles 

The  LIF  signal  in  the  linear  regime  is  proportional  to  both 
the  probed  state  density  and  the  intensity  of  the  probe  beam. 
Therefore,  it  is  at  least  theoretically  possible  to  determine  rela¬ 
tive  ionic  excited  state  number  densities  within  the  plume.  Fig¬ 
ure  1 6  shows  the  LIF  profiles  across  the  exit  plane  with  relative 
signal  intensity  information  retained.  This  shows  that  in  the 
main  peak  of  the  distribution,  the  xenon  excited  state  has  twice 
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Y=8;  Z=-7  mm 
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Fig.  15.  Axial  and  radial  velocity  distributions  (not  deconvolved)  along  the 
nose  of  the  BHT-200-X  for  the  conditions  in  Table  1.  Y  andZ  coordinates  as 
well  as  mean  velocities  are  provided  for  each  data  set.  Note  that  the  VDF 
data  is  area  normalized  and  plotted  with  a  common  scale. 
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the  number  density  at  the  center  than  at  the  channel  edges.  The 
plasma  nearest  the  inner  channel  wall  (Y  =  8  mm)  exhibits  a 
much  larger  proportion  of  low  velocity  ions  than  the  flow  near¬ 
est  the  outer  channel  wall  (Y  =  16  mm).  This  pooling  of  low 
velocity  ions  near  the  inner  channel  wall  may  be  due  to  this 
small  thruster’s  divergent  radial  magnetic  field  which  results  in 
ions  striking  the  inner  walls,  recombining,  and  subsequently 
ionized  to  produce  the  low  velocity  ions,  especially  those  seen 
between  0-5  km/s. 

However,  caution  must  be  taken  in  the  analysis  of  these 
data  in  Fig.  16.  It  is  reasonable  to  state  that  the  ions  may  be 
representative  of  the  local  ions  at  any  single  location,  the 
method  by  which  this  state  is  populated  may  vary  across  the 
exit  plane,  and  certainly  beyond.  The  population  of  this  state  is 
maintained  by  the  simultaneous  processes  of  radiative  decay 
from/to  higher/lower  states,  and  collisional  excitation/de- 
excitation  via  electron  or  heavy  particle  (ion  or  neutral)  colli¬ 
sions.  Therefore  caution  must  be  taken  in  this  sort  of  analysis  in 
general,  but  specifically  here  where  the  plasma  densities  are  so 
variable. 


AIAA-2006-4991 

plume  region  as  the  ions  travel  downstream.  The  axial  VDF 
narrows  due  to  the  external  addition  of  approximately  50  eV  to 
the  flow.  The  radial  VDF  narrows  due  to  the  divergence  of  the 
flow.  Beyond  the  nose  cone,  the  axial  VDF  develops  a  low  en¬ 
ergy  population.  In  a  similar  fashion,  the  radial  VDF  develops  a 
high  velocity  tail.  Both  these  effects  are  most  likely  due  to  cross 
annular  flow,  momentum  exchange  collisions,  and  the  resultant 
mixing  of  velocity  distributions.  Charge  exchange  is  being  in¬ 
vestigated  as  possible  source  of  the  slow  moving  ions,  and  there 
appears  to  be  evidence  for  charge  exchange  collisions  in  the 
near  plume  or  within  the  thruster. 

The  axial  VDF  appears  to  evolve  differently  near  the  inner 
insulator  wall  than  at  the  channel  center.  The  region  near  the 
inner  wall  has  a  higher  proportion  of  low  velocity  ions  in  the 
VDF  and  is  accelerated  25  eV  less  than  the  centerline  case.  A 
large  fraction  of  the  ions  appear  to  be  created  near  the  exit 
plane  at  the  inner  channel  insulator  wall.  The  source  of  the  ions 
is  not  known,  but  from  the  radial  VDF  data,  it  is  suspected  that 
the  slow  ions  may  be  the  result  of  ion  recombination  at  the  in¬ 
ner  wall,  followed  by  local  ionization. 


Fig.  16.  Exit  plane  (1  =  -7  mm)  VDF  contour  and 
surface  plot  scaled  to  LIF  signal  strength  (5d[4ha 
xenon  ion  excited  state).  Note  the  region  of  high 
density  of  low  velocity  ions  along  the  inner  channel 
(Y  =  8  mm). 


Conclusions 

This  work  has  presented  a  preliminary  analysis  of  the  LIF 
derived  velocity  distribution  functions  (VDF)  of  the  5d[4]7/2 
excited  state  of  xenon  in  a  200  W  BHT-200-X3  laboratory  Hall 
thruster.  The  ionic  VDF  evolves  considerably  within  the  near 


Finally,  a  crude  estimate  of  the  relative  densities  at  the  exit 
plane  is  made.  It  appears  that  there  are  a  factor  of  4  more  fast 
moving  ions  in  the  center  of  the  acceleration  channel  than  at  the 
inner  and  outer  regions.  The  relative  width  of  the  VDF  axial 
main  peak  is  near  constant  as  is  the  peak  velocity.  This  indi¬ 
cates  that  these  ions  are  all  created  in  the  same  potential  line 
within  the  thruster.  The  distribution  of  low  velocity  ions  within 
the  VDF  profiles  is  much  stronger  in  the  inner  half  of  the  annu¬ 
lar  acceleration  channel,  possibly  related  to  the  radial  magnetic 
field  divergence  of  this  relatively  small  Hall  thruster. 

Future  Work 

These  results  show  that  the  velocities  of  the  ions  acceler¬ 
ated  by  a  Hall  thruster  can  not  be  well  characterized  by  a  single 
value.  We  believe  that  the  direct  comparison  of  VDF  data  will 
provide  the  best  method  to  validate  numerical  methods.  We 
hope  to  compare  these  data  with  those  derived  from  various 
numerical  models.  Although  bulk  velocities  have  been  com¬ 
pared  in  the  past  [20],  we  are  only  aware  of  one  direct  compari¬ 
son  between  of  VDF  data  derived  from  experimental  and  nu¬ 
merical  means  [21].  The  results  of  this  study  show  that  small 
changes  in  model  parameters  can  produce  dramatic  changes  in 
the  predicted  plasmas.  Although  the  models  appear  able  to 
match  portions  of  the  flow  field,  it  does  not  appear  possible  to 
predict  a  Hall  thruster  ion  velocity  field  a  priori.  Significant 
experimental  inputs  are  required. 

It  would  be  advantageous  to  continue  this  analysis  to  the 
interior  of  the  Hall  thruster.  Many  Hall  thruster  models  use  ex¬ 
perimental  current-voltage  characteristics  to  approximate  the 
electron  conduction  through  radial  magnetic  field.  Therefore,  it 
will  be  possible  to  refine  these  models  with  VDF  data.  Due  to 
the  geometry  of  the  BHT-200-X3  Hall  thruster  it  is  possible  to 
cut  a  slot  to  provide  internal  optical  access  for  the  acquisition  of 
axial  VDF  data  using  the  current  test  apparatus. 

It  is  also  important  to  continue  the  development  of  the  de- 
convolution  routines  into  our  data  reduction  methodology.  We 
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hope  to  develop  a  first  order  Tikhonov  regularization  as  well  as 
tune  the  data  collection  methods  to  increase  our  SNR  to  pro¬ 
duce  clearer  data  for  analysis. 
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